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Abstract 

We discuss whether the NuTeV anomaly can be explained, compatibly with all other 
data, by QCD effects (maybe, if the strange sea is asymmetric, or there is a tiny 
violation of isospin), new physics in propagators or couplings of the vector bosons 
(not really), loops of supersymmetric particles (no), dimension six operators (yes, for 
one specific SU(2)L-invariant operator), leptoquarks (not in a minimal way), extra 
U(l) gauge bosons (maybe: an unmixed Z' coupled to B — also increases the 
muon g — 2hY about 10“® and gives a ‘burst’ to cosmic rays above the GZK cutoff). 


1 Introduction 


The NuTeV collaboration has recently reported a ~ 3 ct anomaly in the NC/CC ratio of deep- 
inelastic t'^-nucleon scattering. The effective coupling to left-handed quarks is found to be about 
1% lower than the best fit SM prediction. 

As in the case of other apparent anomalies (e.g. e'/e ii , the muon g — 2 |, |, |], atomic 
parity violation ^], and another puzzling NuTeV result concerning dimuon events to cite only 
the most recent cases) one should first worry about theoretical uncertainties, mainly due to QCD, 
before speculating on possible new physics. After reviewing in section 2 the SM prediction for the 
NuTeV observables, in section 3 we look for SM effects and/or uncertainties which could alleviate the 
anomaly. In particular, we investigate the possible effect of next-to-leading order QCD corrections 
and consider the uncertainties related to parton distribution functions (PDFs). We notice that a small 
asymmetry between strange and antistrange in the quark sea of the nucleon, suggested by z/Af deep 
inelastic data [10|, could be responsible for a significant fraction of the observed anomaly. We also 
study the effect a very small violation of isospin symmetry can have on the NuTeV result. 


^On leave from INFN, Sezione di Torino, Italy 

^On leave from dipartimento di Fisica deH’Universita di Pisa and INFN. 


1 




SM 

fermion 

U(l)v 

SU(2)i 

SU(3)c 

U'^ 


2 

3 

1 

3 


= d% 

1 

3 

1 

3 



1 

1 

1 

L = 

= (z^L,ei) 

1 

2 

2 

1 

Q-- 

= {uL,dL) 

1 

6 

2 

3 


Z couplings 

9L 

9R 


1 

2 

0 

e,T,T 

-i 

•sw 

u, c, t 

1 2 „2 

2 „2 

2 3'®W 

“3^W 

d, s, b 

_i _L lo2 

2 ^ 3'^W 

1.2 

3'®W 


Table 1: The SM fermions and their Z couplings. 


Having looked at the possible SM explanations, and keeping in mind that large statistical fluc¬ 
tuations cannot be excluded, we then speculate on the sort of physics beyond the SM that could be 
responsible for the NuTeV anomaly. We make a broad review of the main mechanisms through which 
new physics may affect the quantities measured at NuTeV and test them quantitatively, taking into 
account all the constraints coming from other data. 

We take the point of view that interesting models should be able to explain a significant fraction 
of the anomaly. According to this criterion, we consider new physics that only affects the propagators 
(section ^ or gauge interactions (section |5[) of the SM vector bosons, looking at the constraints 
imposed on them by a global fit to the electroweak precision observables. Many models can generate a 
small fraction of the observed discrepancy (see e.g. |^), but it is more difficult to explain a significant 
fraction of the anomaly. In section ^ we consider the case of the minimal supersymmetric SM (MSSM) 
and look at possible MSSM quantum effects. In section we turn to lepton-lepton-quark-quark 
effective vertices, focusing on the most generic set of dimension 6 operators. We find that very few 
of them can fit the NuTeV anomaly (in particular, only one SU(2)i-invariant operator). In section]^ 
and ^ we study how these dimension six operators could be generated by exchange of leptoquarks or 
of extra U(l) gauge bosons. Finally, we summarize our findings in section 10. 


2 The SM prediction 

Tree level 

In order to establish the notation and to present the physics in a simple approximation, it is useful 
to recall the tree-level SM prediction for neutrino-nucleon deep inelastic scattering. The zz^-quark 
effective Lagrangian predicted by the SM at tree level is 

.ifeflF = -2^/2GF{[iy^J,'yahL][dL'y°‘UL] h.C.) - 2V2Gf^ gAq[i^^i'lal^^i][qA'y'^qA] 

A,q 

where A = {L, R}, q = {u, d,s,...} and the Z couplings gAq are given in table ^ in terms of the weak 
mixing angle sw = sin^w 

It is convenient to define the ratios of neutral-current (NC) to charged-current (CC) deep-inelastic 
neutrino-nucleon scattering total cross-sections R^, Rp. Including only first generation quarks, for an 
isoscalar target, and to leading order, these are given by 

^ _ ainM^nX) _{3gl+gl)q+{3gl+gl)q _ ^ 

^ ainM - /rV) “ ^q -^L + rgR (1) 

^ _ a{uM^uX) {3gl + gl)q+i3gl+gl)q ^ , 1^2 

-Kt/ = . r ----WZ- - 9 l + - 9 r , [ 2 ) 

a{uX ^ fiX) q + 3q r 
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where q and q denote the second moments of quark or antiquark distributions and correspond to the 
fraction of the nucleon momentum carried by quarks and antiquarks, respectively. For an isoscalar 
target, q = {u + d)/2, and we have defined 


a{vN fJ.X) _ 3q + q 
a{vM —> fiX) 3q + q 


and 


ai = aiu +9Ld = \- sin^ ^ sin^ %, 


9R^9Ru + 9Rd = gSin %. 


(4) 


The observables and R^^ measured at NuTeV differ from the expressions given in eq. i). 
On the theoretical side this is due to contributions from second-generation quarks, and because of 
QCD and electroweak corrections. On the experimental side, this is because total cross-sections can 
only be determined up to experimental cuts and uncertainties, such as those related to the spectrum 
of the neutrino beam, the contamination of the beam by electron neutrinos, and the efficiency of 
NC/CC discrimination. Once all these effects are taken into account, the NuTeV data can be viewed 
as a measurement of the ratios between the CC and the NO squared neutrino effective couplings. The 
values quoted in Q are 


gl = 0.3005 ± 0.0014 and g% = 0.0310 ± 0.0011, 


(5) 


where errors include both statistical and systematic uncertainties. 

The difference of the effective couplings g\ — g\ (‘Paschos-Wolfenstein ratio’ ||^) is subject to 
smaller theoretical and systematic uncertainties than the individual couplings. Indeed, using eq. i) 
we get 


-Rpw 


Ru ^ Rv 
1 — r 


—> vX) — a{vM vX) ^ 2 

ex) - a{uM ex) 


- - sin2(9w, 


( 6 ) 


which is seen to be independent of q and q, and therefore of the information on the partonic structure 
of the nucleon. Also, i2pw is expected to be less sensitive to the various corrections discussed above. 


Electroweak corrections and the SM fit 


The tree level SM predictions for gi and gR get modified by electroweak radiative corrections. These 
corrections depend on the precise definition of the weak mixing angle, and we therefore adopt the 
on-shell definition |T5| and define = 1 — One then obtains the following expressions for 

aUm 
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z '' 2 A 1 z 

9r — gP , 


(7) 


where, also including the most important QCD and electroweak higher order effects []T| 


p ^ 1.0086 + 0.0001(Mt/GeV- 175)- 0.0006ln(m/,/100GeV), 
k « 1.0349 + 0.0004(Mt/GeV - 175) - 0.0029ln(m,j/100 GeV) 


for values of the top mass, Mt, and of Higgs mass, nih, not far from 175 and 100 GeV. Additional 
very small non-factorizable terms are induced by electroweak box diagrams.|^ 

^QED radiative corrections depend sensitively on the experimental setup and are taken into account in the NuTeV 
analysis. The charged current amplitude receives a small residual renormalization which depends on the implementation 
of QED corrections and is therefore neglected in our analysis. 
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Figure 1: The SM prediction for (g\-,gji) at 
68,99% CL and the NuTeV determination, 
at 68,90,99% CL. The crosses show how the 
NuTeV central value moves along the PW line 
using different sets of parton distribution fune- 
tions that assume s = s. If s > s the PW line 
is shifted towards the SM prediction. 



Figure 2: The NuTeV measurement is eompared 
with the 68,99% CL ranges of (g\,g\) allowed 
by precision data in classes of extensions of the 
SM that allow generic oblique correetions (al¬ 
most horizontal small green ellipses) or flavour- 
universal and SV{2)L-invariant eorrections of 
the Z eouplings (vertieal medium red ellipses). 


We stress that the SM value of sw depends on Mt and mh, unless only the direct measurements 
of Mw and Mz are used to compute it. In particular, for fixed values of Mt and the W mass can 
be very precisely determined from Gp, a{Mz), and Mz. To very good approximation one has (see 

e.g. El) 


Mw = 80.387 - 0.058ln(m/i/100GeV) - 0.008 ln2(m/,/100 GeV) + 0.0062(M4/GeV - 175). ( 8 ) 


Without including the NuTeV results, the latest SM global fit of precision observables gives Mt = 
176.1 ± 4.3GeV, = 87 II 4 GeV, from which one obtains M]y = 80.400 ± O.OlOGeV []T|, and 

therefore = 0.2226 ± 0.0004. The values of g^ and g\ corresponding to the best fit are 0.3042 and 
0.0301, respectively. The small red ellipses in fig. E| show the SM predictions for g\ and g^ at 68 % 
and 99% GL, while the bigger yellow ellipses are the NuTeV data, at 68 % 90% and 99% GL. While 
gj^ is in agreement with the SM, g'j^ shows a discrepancy of about 2.5a. 


Here we have adopted the on-shell definition of sw because it is well known that with this choice 
the electroweak radiative corrections to g^ cancel to a large extent. In fact, at first order in 6p = p — 1 
and 5k = k — 1, g\ gets shifted by Sgf^ ^ {2 5p — 0.551 5k)g\ and the leading quadratic dependence on 
Mt is the same for 5p and 5ksf^/cf^. Therefore, the top mass sensitivity of g\ is very limited when 
this effective coupling is expressed in terms of = 1 — M^/M^. As leading higher order electroweak 
corrections are usually related to the high value of the top mass higher order corrections cannot 
have any relevant impact on the discrepancy between the SM and NuTeV. 

Within the SM, one can extract a value of from the NuTeV data. This is performed by the 
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NuTeV collaboration using a fit to and the effective charm mass that is used to describe the charm 
threshold. This fit is different from the one that gives The result is = 1.32 ± 0.11 GeV and 

= 0.2276 ± 0.0013 (stat.) ± 0.0006 (syst.) ± 0.0006 (th.) (9) 

-0.00003(Mt/GeV - 175) + 0.00032ln(m;,/100 GeV). 


Here the systematic uncertainty includes all the sources of experimental systematics, such as those 
mentioned above, and it is estimated by means of a Monte Carlo simulation of the experiment. The 
theoretical uncertainty is almost entirely given by QCD corrections, to be discussed in section 3. The 
total uncertainty above is about 3/4 of that in the preliminary result from NuTeV = 0.2255 ± 
0.0021 |]^| and about 1/2 of that in the CCFR result = 0.2236 ± 0.0035 |]T|. 

Alternatively, by equating eq. (|^ to the NuTeV results for we hnd 


s^(NuTeV) = 0.2272 ± 0.0017 ± 0.0001 (top) ± 0.0002 (Higgs). 


( 10 ) 


The central value and the errors have been computed using the best global fit for Mt and mh- Eq. (p^) 
has a slightly larger error but is very close to eq. (|9|). An additional difference between the two 
determinations is that eq. is based on a up-to-date treatment of higher order effects. Notice that 
the NuTeV error is much larger than in the global fit given above, from which eq. (P) differs by about 
3(7 and eq. (|I^) by about 2.Ga. The NuTeV result for can also be re-expressed in terms of Mw- 
If we then compare with Mw = 80.451 ± 0.033 GeV from direct measurements at LEP and Tevatron, 
the discrepancy is even higher: more than 3cr in both cases. The inclusion of the NuTeV data in a 
global fit shifts the preferred nih value very slightly, but worsens the fit significantly {x^ = 30 for 14 
degrees of freedom). 

Even without including NuTeV data, the global SM ht has a somewhat low goodness-of-ht, 8% if 
naively estimated with a global Pearson test. The quality of the fit becomes considerably worse if 


only the most precise data are retained |20|. Indeed, among the most precise observables, the leptonic 
asymmetries measured at LEP and SLD and Mw point to an extremely light Higgs, well below the 
direct exclusion bound > 115 GeV, while the forward-backward hadronic asymmetries measured at 
LEP prefer a very heavy Higgs (for a detailed discussion, see |^, ^]). The effective leptonic couplings 
measured by the hadronic asymmetries differ by more then 3 (t from those measured by purely leptonic 
asymmetries. Therefore, the discrepancy between NuTeV and the other data depends also on how 
this other discrepancy is treated. For instance, a fit which excludes the hadronic asymmetries has a 
satisfactory goodness-of-ht, but ruh = 40 GeV as best ht value. In this case, the SM central values for 
9l r 0.3046 and 0.0299, and differ even more from the NuTeV measurements. On the other hand, 
even a very heavy Higgs would not resolve the anomaly: to explain completely the NuTeV result rrih 
should be as heavy as 3 TeV, deep in the non-perturbative regime. The preference of the NuTeV 
result for a heavy Higgs is illustrated in hg. where we display the point corresponding to the SM 
predictions with mh = 500 GeV and mt = 175 GeV. This is suggestive that, as will be more clearly 
seen in the following, the central value of NuTeV cannot be explained by radiative corrections. 


3 QCD corrections 

Most of the quoted theoretical error on the NuTeV determination of is due to QCD effects. 
Yet, this uncertainty does not include some of the assumptions on which the Paschos-Wolfenstein 
relation, eq. (^), is based. Hence, one may ask: first, whether some source of violation of the Paschos- 
Wolfenstein relation which has not been included in the experimental analysis can explain the observed 
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discrepancy, and second, whether some of the theoretical uncertainties might actually be larger than 
estimated in [||. 

A full next-to-leading order (NLO) treatment of neutrino deep-inelastic scattering is possible. 


since all the relevant coefficient functions are long known |^|. If no assumption on the parton content 
of the target is made, including NLO corrections, the Paschos-Wolfenstein ratio eq. @ becomes 


o „2 ,.2 -d ) + {c - s ) 

ripw — 9 l ~ 9 r ^ -qZ- 


“ 9ru) + 2^9Ld ~ 9Rd) 


+ 


+^i9l - 9A{IsA - 5C3)| + 0(Q-)-2 


( 11 ) 


The various quantities which enter eq. ([TI|) are defined as follows: a* is the strong coupling; dC^ = 
6C^ = C^- C* is the the second moment of the next-to-leading contributions to the 
quark coefficient functions for structure function F*; q~ = q — q; Q~ = {u~ + d~)/2; u, d, and so on 
are second moments of the corresponding quark and antiquark distributions. We have expanded the 
result in powers of 1 /Q~, since we are interested in the case of targets where the dominant parton 
is the isoscalar Q~ = {u~ + d~)f2. Equation ( 0 ) shows the well-known fact that the Paschos- 
Wolfenstein relation is corrected if either the target has an isotriplet component (i.e. u ^ d) or sea 
quark contributions have a C-odd component (i.e. 7 ^ 0 or c“ 7 ^ 0). Furthermore, NLO corrections 

only affect these isotriplet or C-odd terms. 

Let us now consider these corrections in turn. Momentum fractions are scale dependent; in the 
energy range of the NuTeV experiment Q~ « 0.18 [^], with better than 10% accuracy, so that 
^idLu ~ 9 ru ) + h^9Ld ~ 9Rd) /2” ~ l'^' Hence, a value {u~ —d~) + {c~ — s~) ^ —0.0038 is required 
to shift the value of by an amount equal to the difference between the NuTeV value central eq. 
and the global SM fit. 

The NuTeV experiment uses an iron target, which has an excess of neutrons over protons of about 
6 %. This violation of isoscalarity is however known to good accuracy, it is included in the NuTeV 
analysis, and it gives a negligible contribution to the overall error. A further violation of isoscalarity 
could be due to the fact that isospin symmetry is violated by the parton distributions of the nucleon, 
i.e. 7 ^ and u'^ 7 ^ d^. This effect is considered by NuTeV Q, but not included in their analysis. 
Indeed, isospin in QCD is only violated by terms of order {niu — md)/^, and thus isospin violation 
is expected to be smaller than 1% or so at the NuTeV scale (where the scale dependence is rather 
weak) |24|. However, if one were to conservatively estimate the associated uncertainty by assuming 
isospin violation of the valence distribution to be at most 1 % (i.e. {u~ — d~)/Q~ < 0 . 01 ), this would 
lead to a theoretical uncertainty on of order As^ = 0.002. This is a more than threefold increase 
in theoretical uncertainty, which would rather reduce the significance of the NuTeV anomaly. 

A C-odd second moment of heavier flavours, 7 ^ 0 or c“ 7 ^ 0 is not forbidden by any symmetry 
of QCD, which only imposes that the first moments of all heavy flavours must vanish. Neither of 
these effects has been considered by NuTeV. A nonzero value of c~ appears very unlikely since the 
perturbatively generated charm component has c“ = 0 identically for all moments, and even assuming 
that there is an ‘intrinsic’ charm component (i.e. c 7 ^ 0 below charm threshold due to nonperturbative 
effects) it is expected to have vanishing c~ ||^ for all moments. On the contrary, because the relevant 
threshold is in the nonperturbative region, the strange component is determined by infrared dynamics 
and there is no reason why s~ = 0. In fact, explicit model calculations [26| suggest s~ ~ 0.002. 
Whereas such a C-odd strange component was at first ruled out by CCER dimuon data |^]^, a 


^Deep-inelastic charm production events are know as dimuon events because their experimental signature is a pair 
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subsequent global fit to all available neutrino data found evidence in favor of a strange component 
of this magnitude and sign |^], and showed that it does not necessarily contradict the direct CCFR 
measurement. More recent measurements [^] confirm the CCFR results in a wider kinematic region; 
however, the quantitative impact of these data on a global fit is unclear. Even though it is not 
included in current parton sets, a small asymmetry s~ ~ 0.002 seems compatible with all the present 

. Assuming s~ 


experimental information |^]. Assuming s ^ 0.002 as suggested by []R|, the value of measured 
by NuTeV is lowered by about = 0.0026. The corresponding shift of the PW line is displayed in 
fig-in This reduces the discrepancy between NuTeV and the SM to the level of about one and a half 
standard deviations (taking the NuTeV error at face value), thus eliminating the anomaly. 


Since NLO corrections in eq. ( [Tl|) only affect the C-odd or isospin-odd terms, they are in practice 
a sub-subleading effect. Numerically, —45(7^ = 16/9 so NLO effects will merely correct a possible 
isotriplet or C-odd contribution by making it larger by a few percent. Therefore, a purely leading-order 
analysis of Rpw is entirely adequate, and neglect of NLO corrections should not contribute significantly 
to either the central value of extracted from i?pW) nor to the error on it. It is important to 
realize, however, that this is not the case when considering the individual ratios and Up. Indeed, 
NLO corrections affect the leading-order expressions by terms proportional to the dominant quark 
component Q~, and also by terms proportional to the gluon distribution, which carries about 50% of 
the nucleon’s momentum. Therefore, one expects NLO corrections to and Rp to be of the same 
size as NLO corrections to typical observables at this scale, i.e. around 10%. The impact of this on 
the values of and however, is difficult to assess: the NuTeV analysis makes use of a parton 
set which has been self-consistently determined fitting leading-order expressions to neutrino data, so 
part of the NLO correction is in effect included in the parton distributions. A reliable determination 
of g\ and g\ could only be obtained if the whole NuTeV analysis were consistently upgraded to NLO. 
As things stand, one should be aware that the NuTeV determination of g\ and g\, eq. (|5|), is affected 
by a theoretical uncertainty related to NLO which has not been taken into account and which may 
well be non-negligible. This uncertainty is however correlated between g‘^ and g"^, and it cancels when 
evaluating the difference g\ — g\. 

On top of explicit violations of the Paschos-Wolfenstein relation, other sources of uncertainty are 
due to the fact that the experiment of course does not quite measure total cross-sections. Therefore, 
some of the dependence on the structure of the nucleon which cancels in ideal observables such as 
Ri, or Rpw remains in actual experimental observables. In order to estimate these uncertainties, we 
have developed a simple Monte Carlo which simulates the NuTeV experimental set-up. The Monte 
Carlo calculates integrated cross sections with cuts typical of a I'M experiment, by using leading- 
order expressions. Because the Monte Carlo is not fitted self-consistently to the experimental raw 
data, it is unlikely to give an accurate description of actual data. However, it can be used to assess 
the uncertainties involved in various aspects of the analysis. 

We have therefore studied the variation of the result for i?pw as several theoretical assumptions 
are varied, none of which affects the ideal observable Rpw but all of which affect the experimental 
results. First, we have considered the dependence on parton distributions. Although the error on 
parton distributions cannot really be assessed at present, it is unlikely to be much larger than the 
difference between leading-order and NLO parton sets. We can study this variation by comparing the 
CTEQL and CTEQM parton sets |Q. We also compare results to those of the MRST99 set |^. We 
find extremely small variation for i?pw and small variations even for the extraction of Specific 

uncertainties which may affect significantly neutrino cross sections are the relative size of up and down 


of opposite-sign muons. Since they mainly proceed through scattering of the neutrino off a strange quark, they are a 
sensitive probe of the strange distribution. 
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distributions at large x |^, and the size of the strange and charm component |^|. Both have 


been explored by MRST [23|, which produce parton sets where all these features are varied in turn. 
Again, using these parton sets in turn, we find no significant variation of the predicted Ry, Rp, and 
of the extracted If, on the contrary, we relax the assumption s = s, which is implicit in all these 
parton sets, we find a shift of i?pw in good agreement with eq. (D)- This conclusion appears to be 
robust, and only weakly affected by the choice of parton distributions and by the specific x-dependence 
of the s — s difference, provided the second moment of s — s is kept fixed. The lower cut (20GeV) 
imposed by NuTeV on the energy deposited in the calorimeter tends to decrease the sensitivity to the 
asymmetry s~, as it mostly eliminates high x events. However, this effect is relevant only for lower 
energy neutrinos, below about 100 GeV, and should be small in the case of NuTeV. 

The dependence on the choice of parton distributions is shown in fig. || where blue x (red +) 
crosses correspond to MRST99 (CTEQ) points. We cannot show a NuTeV value, because we could 
not access the parton set used by NuTeV. The results are seen to spread along the expected PW 
line. The intercept of this line turns out to be determined by the input value of g^ — 9%, and to 
be completely insensitive to details of parton distributions. However, it should be kept in mind that 
inclusion of NLO corrections might alter significantly these results, by increasing the spread especially 
in the direction along the PW line, for the reasons discussed above. 

Finally, we have tried to vary the charm mass, and to switch on some higher twist effects (specifi¬ 
cally those related to the nucleon mass). In both cases the contributions to the uncertainty which we 
find are in agreement with those of NuTeV. 


4 Oblique corrections 


After our review of the SM analysis, let us proceed with a discussion of possible effects of physics 
beyond the SM. We first concentrate on new physics which is characterized by a high mass scale and 
couples only or predominantly to the vector bosons. In this case its contributions can be param¬ 
eterized in a model independent way by three (oblique) parameters. Among the several equivalent 

Many models of physics beyond the SM can be studied 


parameterizations ||^, we adopt ei,€ 2 ,e 3 
at least approximately in this simple way. 

Generic contributions to ei,e 2 ,e 3 shift g^ according to the approximate expressions 


M 

9l 


= 2.8(5ei - 1.1(563; 


^9r 

9r 


= -0.9(5ei -F 3.7(563. 


( 12 ) 


Of course, the 6 j parameters are strongly constrained by electroweak precision tests. In order to see 
if this generic class of new physics can give rise to the NuTeV anomaly, we extract the 6 * parameters 
directly from a fit to the electroweak data, without using the SM predictions for them. We use the 
most recent set of electroweak observables, summarized in table |^, properly taking into account the 
uncertainties on aem(Mz) and as{Mz)- The result is a fit to the 6 * very close to the one reported 
in 1^ which we use in eq. (12) after normalizing to the SM prediction at a reference value of m/j. 
The ellipses corresponding to the 68 % and 99% CL are displayed in fig. Q (green, almost horizontal 
ellipses). They are centred roughly around the SM best fit point, because the SM predictions for 
6 i and 63 for rrih ~ 100 GeV are in reasonable agreement with the data (see also section 2). The 
difference between the best ht point and the light Higgs SM prediction for (fl'£,fl'^) is much smaller 
than the NuTeV accuracy. Notice that, as mentioned in section 2, excluding the hadronic asymmetries 
from the fit would make an oblique explanation even harder. 

Our conclusion is that oblique corrections cannot account for the NuTeV anomaly, as they can 
only absorb about ~ Icr of the ~ 3a NuTeV anomaly. 







Gf 

= 

1.16637 lO-VGeV^ 


Mz 

= 

91.1875 GeV 


mt 

= 

(174.3 ± 5.1) GeV 


as{Mz) 

= 

0.119 ±0.003 


aVrtiMz) 

= 

128.936 ± 0.046 


M\y 

= 

(80.451 ± 0.033) GeV 

1 .8-cj 

Tz 

= 

(2.4952 ± 0.0023) GeV 

—0.4-cj 

(Th 

= 

(41.540 ± 0.037)nb 

1 .6-cj 

Ri 

= 

20.767 ± 0.025 

1 .1-CJ 

Rb 

= 

0.21646 ± 0.00065 

1 .1-CJ 

Rc 

= 

0.1719 ±0.0031 

-0.1-cj 

^LR 

= 

0.1513 ±0.0021 

1 .6-cj 


= 

0.922 ± 0.02 

—0.6-cj 

^LR 

= 

0.670 ± 0.026 

0 .1-cj 

^FB 

= 

0.01714 ±0.00095 

0 .8-cj 

^FB 

= 

0.099 ± 0.0017 

-2.8-cj 

/ic 

^FB 

= 

0.0685 ± 0.0034 

-1.7-cj 

Qw 

= 

-72.5 ±0.7 

0 .6-cj 


Fermi constant for // decay 

pole Z mass 

pole top masss 

strong coupling 

electromagnetic coupling 

pole W mass 

total Z width 

ee hadronic cross section at Z peak 
Viz ^ hadrons)/r(Z ^ 

V{Z bb)/V{Z hadrons) 

V{Z —> cc)/T{Z —> hadrons) 

Left/Right asymmetry in ee 
LR Forward/Backward asymmetry in ee —> bb 
LR FB asymmetry in ee —> cc 
Forward/Backward asymmetry in ee ^ ll 
Forward/Backward asymmetry in ee ^ bb 
Forward/Backward asymmetry in ee ^ cc 
atomic parity violation in Cs 


Table 2: The electroweak data included in our fit J^/ . The second column indicates the discrepancy 
with respect to the best SM fit. 


5 Corrections to gauge boson interactions 

We now discuss whether the NuTeV anomaly could be explained by modifying the couplings of the 
vector bosons. This possibility could work if new physics only affects the hZi' couplings, reducing 
the squared Vf^Zu^ coupling by (1.16 ± 0.42)% |]|]. This shift is consistent with precision LEP data, 
that could not measure the uZu couplings as accurately as other couplings (no knowledge of the 
LEP luminosity is needed to test charged lepton and quark couplings), and found a Z ^ i/h rate 
(0.53 ± 0.28)% lower than the best-fit SM prediction. We could not construct a model that naturally 
realizes this intriguing possibility, because precision data test the flZfj. and couplings at the per 

mille accuracy. This generic problem is best understood by considering explicit examples. We first 
show that models where neutrinos mix with some extra fermion (thereby shifting not only the z/Zzz 
coupling, but also the iWi' coupling) do not explain the NuTeV anomaly. Next, we discuss why a 
model where the Z mixes with some extra vector boson (thereby shifting not only the hZi' coupling, 
but also the Z couplings of other fermions) does not explain the NuTeV anomaly. 


Models that only affect the neutrino couplings 

This happens e.g. in models where the SM neutrinos mix with right-handed neutrinos (a 1% mixing 
could be naturally obtained in extra dimensional models or in appropriate 4-dimensional models [^]). 
By integrating out the right-handed neutrinos, at tree level one obtains the low energy effective 
lagrangian 

Jf eff = Jf SM + ^ij 2V2GFiH^Li)id{HLj), (13) 
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where i,j = {e,/x,r}, Li are the lepton left-handed doublets, H is the Higgs doublet, and = 
e*ji are dimensionless couplings. This peculiar dimension 6 operator only affects neutrinos. After 
electroweak symmetry breaking, it affects the kinetic term of the neutrinos, that can be recast in the 
standard form with a redefinition of the neutrino field. In this way, the ViZvi and the UiWii couplings 
become respectively 1 — eu and 1 — ejj/2 lower than in the SM {ea is positive: gauge couplings of 
neutrinos get reduced if neutrinos mix with neutral singlets). The NuTeV anomaly would require 
C/i/i = 0.0116 ± 0.0042. However, a reduction of the eWve and flWv^ couplings increases the muon 
lifetime, that agrees at about the per-mille level with the SM prediction obtained from precision 
measurements of the electromagnetic coupling and of the W and Z masses. Assuming that no other 
new physics beyond the extra operator in eq. (13) is present, from a fit of the data in table § we find 
that a flavour-universal eu is constrained to be eu = (0 ± 0.4)10“^. This bound cannot be evaded 
with flavour non universal corrections, that are too strongly constrained by lepton universality tests 

In conclusion, Cnn can possibly generate only a small fraction of the NuTeV 


in T and vr decays 
anomaly. 

In principle, the strong bound from muon decay could be circumvented by mixing the neutrinos 
with extra fermions that have the same W coupling of neutrinos but a different Z coupling. In practice, 
it is not easy to build such models. 


Models that only affect the Z couplings 

Only the Z couplings are modihed, e.g., in models with an extra U(l) Z' gauge boson that mixes with 
the Z boson. The Z' effects can be described by the ZZ' mixing angle, 6, by the Z' boson mass, Mz', 
and by the Z gauge current Jz'- At leading order in small 9 and Mz/Mz', the tree-level low energy 
lagrangian gets modihed in three ways. 

(1) the SM prediction for the Z mass gets replaced by 

(2) the Z current becomes Jz = J^ — OJz'', 

(3) at low energy, there are the four fermion operators generated by Z' exchange, beyond the ones 
generated by the W± and Z bosons: 


■Z’es{E <C Mz, Mz') = 


Jw+Jw- 

Ml 


t2 t2 

I 




As discussed in section (1) cannot explain the ~ 1% NuTeV anomaly. Here we show that the 
same happens also for (2): the Z couplings are constrained by LEP and SLD at the per-mille level, 
and less accurately by atomic parity violation data, as summarized in table However, the less 
accurate of these data have ~ 1% errors, and present some anomalies. The Z ^ vv rate and the 
Forward/Backward asymmetries of the h and c quarks show a few-u discrepancy with the best-ht 
SM prediction. But the Z ^ bb and Z ^ cc branching ratios agree with the SM. The best SM fit, 
including also the NuTeV data Q, has ~ 30 with 14 d.o.f. In this situation, it is interesting to 
study if these anomalies could have a common solution with Z couplings about 1% different from the 
SM predictions. We therefore extract the Z couplings directly from the data, without imposing the 
SM predictions for them. This kind of analysis has a general interest. Since we are here concerned 
with the NuTeV anomaly, we apply our results to compute the range of consistent with the 

electroweak data. We recall that both neutrino and quark couplings enter in determining and gj^. 

We assume that the Z couplings are generation universal and S\J{2)l invariant as in the SM: 
we therefore extract from the data the 5 parameters gg, gjj, go, 9 l and gE that describe the Z 
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couplings to the five kinds of SM fermions listed in table ||. In the context of Z' models, this amounts 
to assume that the Z' has generation-universal couplings that respect 311(2)/,. This assumption of 
SU(2)/, invariance is theoretically well justified, although one could possibly invent some non minimal 
model where it does not hold. On the contrary, the universality assumption only has a pragmatic 
motivation; we cannot make a fit with more parameters than data. 

We obtain the result shown by the large red ellipse on the right side of fig. This generic class 
of models gives a best ht value close to the SM prediction. Although the error is much larger than 
in a pure SM fit, it does not allow to cleanly explain the NuTeV anomaly. We find that the global 
can be decreased by about 4 with respect to a SM fit: taking into account that we have five more 
parameters this is not a statistically significant reduction^ in agreement with old similar analyses 


One could generalize this analysis in several directions. For example, new physics could shift the 
on-shell Z couplings tested at LEP and SLD differently from the low-energy Z couplings relevant for 
NuTeV. Alternatively, there could be flavour dependent shifts of the Z couplings. This happens e.g. 
in the model considered in |n], where it is suggested that the NuTeV anomaly could be reproduced 
by a mixing between the Z boson with a Z' boson coupled to the lepton flavour numbers — Lr- 
However, this mixing also shifts the couplings of charged r and ^ leptons, that are too precisely tested 
by LEP and SLD to allow for a significant fraction of the NuTeV anomaly. 


6 Loop effects in the MSSM 


It is well known that supersymmetric contributions to the electroweak precision observables decouple 
rapidly. Under the present experimental constraints it is very difficult to find regions of parameter 
space where radiative corrections can exceed a few per-mille. Explaining the NuTeV anomaly (a 1.2% 
discrepancy with the SM prediction for g'f^) with low-energy supersymmetry looks hopeless from the 
start. Moreover, the dominant contributions to ei in the MSSM are always positive |^]. It then follows 
from eqs. (12) that, in order to explain at least partially the measured value of the supersymmetric 
contributions to €3 should be positive and of 0(1%). 


An interesting scenario which can be easily investigated is the one recently proposed in [^] to 
improve the global fit to the electroweak data. As the main contributions of squark loops would 
be a positive shift in ei, all squarks can be assumed heavy, with masses of the order of one TeV. 
Relatively large supersymmetric contributions are then provided by light gauginos and sleptons and 
can be parameterized in terms of only four supersymmetric parameters (tan /3, the Higgsino mass //, 
the weak gaugino mass M 2 , and a supersymmetry-breaking mass of left-handed sleptons). The oblique 
approximation used in section Q is not appropriate for light superpartners (sneutrinos can be as light 
as 50GeV). We therefore consider the complete supersymmetric one-loop corrections in this scenario 
(see [^] and refs, therein). Taking into account the various experimental bounds on the chargino and 
slepton masses, we hnd the potential shifts to gj^ ^ shown in fig. They are small and have the wrong 
sign. Low-energy supersymmetric loops cannot generate the NuTeV anomaly. 


^On the contrary, if we allow for generation universal but non SU(2)_L-invariant corrections of the Z couplings to ul 
and di and to v and cl, we get a statistically significant reduction in the global ~ 21 with 7 more parameters), 

due to the fact that various anomalies, including of course the NuTeV anomaly, can be explained in this artificial context. 
Without including the NuTeV data, the best fit regions in the {g\,g^ plane are shifted towards the NuTeV region. 
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Figure 3: Shifts in g\ ^ in the supersymmet¬ 
ric scenario described in the text (light slep- 
tons and gauginos) for points in the parame¬ 
ter space which are not excluded by experimental 
constraints. 


Figure 4: Deviations from the SM prediction 
that can be induced by minimal leptoquarks (con¬ 
tinuous lines) and by an extra B — 2>L^ gauge bo¬ 
son (dashed red line) without conflicting by more 
than la with other bounds. 


7 Non renormalizable operators 

Non renormalizable operators parameterize the effects of any new physics too heavy to be directly 
produced. As discussed in sections ^ and new physics that affects the Z, W± propagators or couplings 
cannot ht the NuTeV anomaly without some conflict with other electroweak data. We now consider 
dimension six lepton-lepton-quark-quark operators that conserve baryon and lepton number. 

We start from a phenomenological perspective, with SU(3)(8>U(l)em invariant four fermion vertices, 
and determine which vertices could explain the NuTeV anomaly without conflicting with other data. 
Then we next consider which SU(2)2, invariant operators generate the desired four-fermion vertices. 
In the next sections, we will discuss new particles whose exchange could generate these operators. 

Taking into account Fierz identities, the most generic Lagrangian that we have to consider can be 
written as 


.if 


eflf 


= .ifsM-2\/2Gj 


4Xqrqs i^a^PAij){nPBqs) + 


P ^iiljqrqs^^'^^^ PA(-j){qrlfiuPAgs) 


(14) 


where = f [7/^,7^], Pr,l = (1 =t 75)/2 are the right- and left-handed projectors, q and i are any 
quark or lepton, A,B = {L, R} and e, 6 and t are dimensionless coefficients. In order to explain the 
NuTeV anomaly, new physics should give a negative contribution to g^. This can be accomplished by 


1. reducing the NC z^^-nucleon cross section. The operators given by new physics must contain 
left-handed first generation quarks. 
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2 . increasing the CC z^^-nucleon cross section. The quarks do not need to be left-handed, and the 
quark in the final state needs not to be of the first generation. 

We now show that ‘scalar’ operators (the ones with coefficients 5) cannot explain NuTeV, left- 
handed ‘vector’ operators (with coefficient can realize the first possibility and ‘tensor’ operators 
(with coefficient t) perhaps the second one. 

The scalar operators with coefficient 5 contribute to the charged current. In order to accom¬ 
modate the NuTeV anomaly, these operators should appear with a relatively large coefficient <5 >0.1, 
since their contribution to CC scattering has only a negligible interference with the dominant SM 
amplitude. The interference is suppressed by fermion masses (for first generation quarks) or by CKM 
mixings (if a quark of higher generation is involved). For first generation quarks, this value of 
is inconsistent with When new physics-SM interference is included in this ratio, it becomes |^| 


^ BR(7r ^ eue) ^ gM 
" - BR(7r ^ 


1 — 9(f-^ - — 1 — 

^Up,ud ^eue,ud) 


2m 


_eiP _ 

+ mrf) 


(15) 


The measured value, Rt^ = (1.230 ± 0.004) x 10 Q, agrees with the SM prediction 

mg(m^ — TTig)^ 


,0 


^SM ^ 


(m^ — 


(1 - 16.2 




TT 


which implies Sf^y^^d ^ 10“^. Furthermore, scalar operators which produce a s,c or 6 quark in the 
final state also cannot explain the NuTeV anomaly. The values of 5 required would be in conflict with 
upper bounds on FCNC meson decays such as —> vr+^r/i, D~^ —> fifi, and —> fifl. 


Vector operators can possibly generate me in urev anumaiy 
first generation quarks, the operators in eq. (14) shift as 










— 2 [gii, 


LL 


"Uu,Vu,UU 


, „ .LL 
uu 1 “ yLd ^p^y^dd 


— 9l ^ 


XL _ n 


(16) 


The CC term, cannot alone fit the NuTeV anomaly without overcontributing to the tt 

decay. In principle, one could allow for cancellations between different contributions to R-j^ in eq. (|l^). 

However LEP (and bounds from atomic parity violation |^^) exclude the simplest possibility, 
LL _ _ LL _ 

^eueud' 

We now assume that these vector operators are generated by new physics heavier than the max¬ 
imal energy of present colliders (about few hundred GeV), and study the bounds from collider data. 
Operators involving second generation leptons are constrained by the Tevatron; LEP and HERA are 
not sensitive to them^ In the case of vector operators, the Tevatron sets a limit ^0.03 

(g = M is slightly more constrained than q = d, because protons contain more u than d). Assuming 
SU(2)x,-invariance (that relates the with i = {//, and q = {u,d}, see below) the Tevatron 
bound is close but consistent with the value suggested by NuTeV, l^p^y^gql ~ 0.01. 


Tensor operators could possibly produce the NuTeV anomaly via mechanism 2, because vr-decays 
give no bound on tpy^y^d (using only the tt momentum it is not possible to write any antisymmetric 
tensor). Tensor operators have not been studied in [^|, but if they are generated by physics at a scale 


®Both at NuTeV and at CCFR these operators may induce deviations from the SM values of qlXR- Within their 
errors, the CCFR data are consistent with both the SM and with NuTeV, as clearly shown in fig. 3 of [Q. Therefore 
CCFR bounds on non renormalizable operators, reported by the PDG cannot conflict with NuTeV. 
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^ Mz, the value of ~ 0.1 necessary to fit the NuTeV anomaly is within (and probably above) 

the sensitivity of present Tevatron data. Furthermore we do not know how new physics (e.g. exchange 
of new scalar or vector particles could generate only tensor operators, without also generating 
the scalar operators that overcontribute to R.,^. We will therefore focus on vector operators. 


We now consider SU(2)x,-invariant operators. We have shown that the NuTeV anomaly could be 
explained by the four-fermion vertex Only two SU(2)i invariant operators 

operators can generate this vertex 

Olq = [L7^L] [Q7^Q] , O'lq = [h^T^L] [Q^^r^Q ]. 

We left implicit the 811(2)/, indices, on which the Pauli matrices act. Other possible 4 fermion 
operators, with different contractions of the SU(2)/, indices, can be rewritten as linear combinations 
of these two operators. 

The NuTeV anomaly can be fit by Olq if it is present in eff as (—0.024 ± 0.009) 2 -v/2GfOlq, as 
discussed above. The operator 


^'lq = “ IJ-Ll,j.lJ-L][uLlfiUL - dL'^fj.dL] + 2[/iL7^i^^]-h 2[i^^'y^fiL][dL7f,UL] 


also can fit the NuTeV anomaly. However its CC part overcontributes to vr —> giving a contribu¬ 
tion to about 10 times larger than what allowed by R-,^, see eq. ([l^ . 

These operators could be induced e.g. by leptoquark or Z' boson exchange, which we study in the 
following two sections. A critical difference between these possibilities is that leptoquarks must be 

|, whereas a neutral Z' boson could also be lighter than about 
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heavier than about 200 GeV 

lOGeV (see section ^). Leptoquarks are charged and coloured particles that would be pair-produced 
at colliders, if kinematically possible. If the NuTeV anomaly is due to leptoquarks, their effects should 
be seen at run II of the Tevatron or at the LHG. If instead the NuTeV anomaly were due to a weakly 
coupled light Z', it will not show up at Tevatron or LHG. 


8 Leptoquarks 


Leptoquarks are scalar or vector bosons with a coupling to leptons and quarks. In this section, we 
consider leptoquarks which induce baryon and lepton number conserving four-fermion vertices. 


The symmetries of the SM allow different types of leptoquarks, which are listed in |51]. There are 
four leptoquarks that couple to QL, so these are candidates to explain the NuTeV anomaly. They are 
the scalar SU(2)/, singlet (5o) and triplet (5f), and the vector 811(2)/, singlet (Vq^) and triplet 
with interaction Lagrangian 


Aso [QL]So + A5 i [Qr'^LjSf + Xvo + Ay^ + h.c. 

We do not speculate on how the above leptoquarks could arise in specific models. 


(17) 


Gonsider first the scalar Sq. The lower bound on leptoquark masses from the Tevatron is 200 GeV K 
^1, therefore at NuTeV leptoquarks are equivalent to effective operators. Tree level exchange of Sq, 
with mass m and coupling X[QiL 2 ]Sq (1 and 2 are generation indices), induces the four-fermion oper¬ 
ator 

^,« = ^(Qli2)(QlL2)' = ^{OlQ - Oig). 
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LQ 

-Sf eff 


gLL 

Uu.^u.dd 

LL 

LL 

^icpiiypiuu 

^LL_ 

flfidd 

LL LL 

Sq 

^(Olq - ote) 

0 .12a 

-ajl 

-a/2 

0 

0 

a/2 




0 

0 

—a 

0 

0 

Si 

iil|(0te + 30ig) 

0.03a 

-ajl 

-a/2 

0 

0 

-a/2 




0 

0 

0 

—a 

0 

Vq 

S(o'LQ + OLe) 

0.09a 

0 

0 

a 

a 

a 




0 

2 a 

0 

0 

0 

Vi 

- 30ig) 

—0.40a 

0 

0 

a 

a 

—a 




2 a 

0 

0 

0 

0 


Table 3: The four leptoquarks that generate the operators suggested by the NuTeV anomaly. Columns 
two and three are the effective Lagrangian and the contribution to g\ (in terms of a = \)(^\/2^/2G■pm?‘) 
for leptoquarks with SU(2)2, degenerate masses. In the last columns, we give the coefficients of the 
individual four fermion operators, separating the triplet members. 


The sign of the operator is fixed, because the coupling constant is squared. We see that Sq can¬ 
not explain the NuTeV anomaly, because it generates Olq with the wrong sign (it gives a positive 
contribution to gf^), and because it also generates the unwanted operator Ofq. 

In the context of supersymmetric models without i?-parity Sq can be identified with a Dg squark 
of generation g and superpotential interaction X' 2 giL 2 DgQi. It is interesting to explore further the 
possible contributions of i?-parity violating squarks at NuTeV. In supersymmetric models, is 
accompanied by a scalar SU(2)i doublet squark (leptoquark), Q. The exchange of Q only modifies 
the right-handed coupling gn, so that it cannot explain NuTeV by itself. Mixing of right- and left- 
handed squarks generates dimension seven operators. This mixing is usually, but not always, negligibly 
small (e.g. one can consider large tan/3, or non minimal models). The relevant AT = 2 four-fermion 
operators are 


^ijk^mnj [dkPLk'illQnPLTm] — ^ijk^mnj[^kPL^i]['H‘nPL(im ( 18 ) 

These operators cannot account for the NuTeV anomaly: they do not interfere with the SM amplitude 
and contribute to both NC and CC, leading to a positive correction to gi. 

In table we list the effective four-fermion operators, and the contribution to g\, of Sq, Si, Vq 
and Vi. In the .5? eff column we have assumed that the members of triplet leptoquarks are degenerate. 
Only the vector SU(2)x, triplet leptoquark gives a negative contribution to g\. In all cases Olq is 
generated together with the unwanted Ofq operator, that overcontributes to the vr —> decay, as 

discussed in the previous section. These features are also shown in fig. where we plot the deviations 
from the SM prediction induced by the Sq, Si, Vq, Vi leptoquarks imposing that they should not 
overcontribute to Rt^ by more than la. 

In the subsequent columns of table ^ we generalize the effective Lagrangian assuming that SU(2)i 
breaking effects split the triplets in a general way. In this situation, the scalar and vector SU(2)i 
triplet leptoquarks can explain the NuTeV anomaly. In the scalar (vector) case, NuTeV can be fit by 
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reducing the mass of the triplet member that induces ^ factor of \/2. From |^ , 

we expect that such split multiplets are consistent with precision electroweak measurements. 

We conclude that the NuTeV anomaly cannot be generated by SU(2)i singlet or doublet lepto- 
quarks, or by triplet leptoquarks with degenerate masses. However, triplet leptoquarks with carefully 
chosen mass splittings between the triplet members can fit the NuTeV data — and this explanation 
should be tested at Run II of the Tevatron or at LHC. 


9 Unmixed extra U(l) Z' boson 

The sign of the dimension 6 lepton/quark operators generated by an extra Z' vector boson depends 
on the lepton and quark charges under the U(l)^ gauge symmetry. Therefore, with generic charges, 
it is possible to generate a correction to neutrino/nucleon scattering with the sign suggested by the 
NuTeV anomaly. In order to focus on theoretically appealing Z' bosons, we require that 


• Quark and lepton mass terms are neutral under the extra U(l)h We make this assumption 
because experimental bounds on flavour and CP-violating processes suggest that we do not have 
a flavour symmetry at the electroweak scale. 

• The Z' couples to leptons of only second generation. Bounds from (mainly) LEP2 and 
older ee colliders would prevent to explain the NuTeV anomaly in presence of couplings to hrst 
generation leptons. We avoid couplings to third generation leptons just for simplicity. 

• The extra U(l)' does not have anomalies that require extra light fermions charged under the SM 
gauge group. 


The only gauge symmetry that satishes these conditions is H — 3L^ (for related work see ||5^ ), where 
B is the baryon number and is the muon number.^ Under these restrictions, the sign of the Z' 
correction to neutrino/nucleon scattering is fixed, and this Z' allows to ht the NuTeV anomaly. In 
fact, the four-fermions operators generated by Z' exchange are 


=_ ^ _ 

^ 2(M|, - t) [ 

where t is momentum transferred: t 
obtained for (see fig. Q) 


- 9L27/.^2 + 

-20GeV2 at NuTeV. The best fit of 


2 


the NuTeV anomaly is 


M|, — t Ri g^/ 3 TeV. 

We now discuss the experimental bounds on such a Z'. 


(19) 


Collider bounds 


The bounds from Tevatron 


m 


a{pp Z'X sXZs = 1.8 TeV)BR(Z' ^ pg) < 40 fb (95%CL) 

®The neutrino masses suggested by oscillation data do not respect L^. If we allow third generation leptons to be 
charged under the extra U(l)' symmetry, we could have a. B — cLfj, — (3 — c)Lt gauge group (c is a constant). However 
even a ± L-r symmetry would not force a successful pattern of neutrino masses and mixings, that rather suggest a 
flavour symmetry containing Le — L/i “ Lt . 
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and LEP @[| 


-6 


BR(Z —> fj,^Z')YiK{Z' —> < few x 10 

imply that Mz' cannot be comparable to Mz- One needs either a light Z', Mz' ^ 10 GeV, or a heavy 
Z', Mz' ^600 GeV [^|. Perturbativity implies Mz' ^5TeV. 


The anomalous magnetic moment of the muon 


The Z' gives a correction to the anomalous magnetic moment of the muon. Assuming Mz' 
get 


a,, = a 


SM 


+ 


274, 


rriT 


47r2 Ml, 


= a 


SM 


8.4 10 


-10 ( 


3TeV y 
^Mz'/gz') 


we 


At the moment, the measured by the g—2 collaboration Q is slightly higher than the SM prediction, 
ayP — = (20 ± 18)10“^°, if one employs = (697 ± 10)10“^*^ for the hadronic polarization 

contribution Q and = (9 ± 2) 10“^*^ for the light-by-light contribution |^]. One could still prefer 
to estimate it as = (0 ± 10) 10~^^, obtaining a larger discrepancy and error. 

Therefore, if Mz' / gz' ~ 3 TeV, as suggested by the NuTeV anomaly in the heavy Z' case, the Z' 
correction to the g — 2 is comparable to the sensitivity of present experiments. If instead Mz' ^ 5 GeV, 
the Z' that fits the NuTeV anomaly gives a larger correction to the muon g — 2, see eq. (^). For 
example, for Mz' ~ 3 GeV one can fit the central value of — a^. On the other hand, a Z' lighter 
than (1 2) GeV cannot explain the NuTeV anomaly without overcontributing to a^. Similar light 

Z' models were proposed |5^ as a possible source of the discrepancy between and previous SM 
computations |Q]. 


Other bounds 

Quantum corrections generate an unwanted kinetic mixing between the Z' and the SM hypercharge 
boson A light Z' needs a small gauge coupling gz', making these quantum effects negligible. 

The Z' could contribute to the decay of qq mesons into fig. This is negligibly small for gz' ~ 
Mz'/3TeV , unless mz' is very close to the meson mass nigg. There are various cc mesons in the 
few GeV mass range, but P^/ is very narrow, so the Z' is only ruled out in narrow windows mz' ~ 

mgg±Tgg Q. 

The Z' that can fit the NuTeV anomaly does not give significant corrections to rare K, D and 
B decays. Let us consider for example the K~^ —> decay. This is a sensitive probe because 

the dominant SM Z penguins that know that SU(2)i is broken (and can therefore generate the 
FCNC vertex m^[sL^^Z^dL\, with GIM cancellations spoiled by the large top mass) are suppressed 
by the small mixing V^*Vtd. On the contrary, penguin loops of quarks that do not know that U(l)' is 
broken (that therefore only generate the q'^lsL'y^Z'^di] operator, where q is the Z' momentum and the 
GIM cancellations is only spoiled by logarithms of quark masses) are suppressed only by the larger 
Cabibbo mixing Vcd- The Z' suggested by NuTeV gives a negligible correction even to this favorable 
K~^ —> decay 

^The total number of measured Z agrees with the SM prediction and there is no peak in the gfl invariant 

mass. In order to extract a precise bound on the Z' mass from these data one should take into account the experimental 
cuts and resolution. 

®A Z' boson with mass Mz' ~ 100 MeV (this case is not motivated by NuTeV) would mediate the resonant decay 
K+ —> Z' -K^vv, producing an excess of monoenergetic in the rest frame, compatibly with the first 
experimental data M- 
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Ref. 1^] claims a statistically significant hint of dimuon events (three events, versus an estimated 
background of 0.07±0.01 events), possibly generated by some neutral long lived particle with produc¬ 
tion cross section 10 / GeV^ of few GeV mass (see also @). The Z' suggested by NuTeV can 
have the right mass and cross section, but is not sufficiently long lived. However, it could partially 
decay in sufficiently long lived neutral fermions. Extra light neutral fermions are required to cancel 
gravitational anomalies and U(l)'^ anomalies. 


Z' burst and the GZK cutoff 


The Z' gives a narrow resonant contribution to vv scattering, which could perhaps generate ultrahigh 
energy cosmic ray events with E ~ 10^^ GeV. The analogous Z resonance [^] has been considered 
as a possible source of the observed events above the Greisen-Zatsepin-Kuzmin (GZK) cutoff [^]. A 
cosmic ray neutrino that scatters with a nonrelativistic cosmic microwave background neutrino would 
encounter the Z' resonance at the energy 


jpZ' res 


2m^ 


10^0 eV 


Mz' 

3 GeV 


^O.OheV 

rriu 


( 20 ) 


where we have used a neutrino mass suggested by atmospheric oscillation data A resonance at 
a Z' mass of few GeV is more suitable than a resonance on the Z, where a larger incident neutrino 
energy >4 X 10^^ eV would be required, even if neutrinos are as heavy as possible, 

my < 1 eV. The Z burst scenario is problematic because it seems difficult to imagine a cosmological 
source that produces enough very energetic neutrinos without producing, at the same time, too many 
photons. 

Although the Z' requires less energetic cosmic neutrinos than the Z, roughly the same total flux [^| 
is required in the two cases, because the Z and the Z' have comparable energy-averaged cross sections. 
In fact the NuTeV and g — 2 data suggest Vz' /Mz' ~ VzjMz and 


cj{vv ^ Z' ^ f) 


^ /)BR(Z' ^ vv)b{s - Ml) 

Mz' 


where / is any final state. In conclusion, a Z' burst could generate the observed cosmic rays above 
the GZK cutoff more easily than the Z burst. 


10 Summary 


We have studied the possible origin of the NuTeV anomaly. Our main results are: 

• QCD effects. Because of the approximate use of the Paschos-Wolfenstein relation, the dis¬ 
crepancy between the NuTeV result and the SM prediction is fairly independent of which set 
of standard parton distribution functions is employed. This is no longer true if one drops some 
of the simplifying assumptions which are usually made in the PDF fits. A small asymmetry 
between the momentum carried by strange and antistrange quarks in the nucleon, suggested by 


a recent analysis of neutrino data |10|, could explain half of the discrepancy between NuTeV 
and the SM. Such an asymmetry has been set to zero in more recent global parton sets, but the 
value suggested in 0 seems compatible with the data used in these hts. It would be desirable 
to update the analysis of |T^ including more recent di-muon data from NuTeV |^, though 
settling the issue is likely to require more detailed information, such as it could be achieved at 
a neutrino factory p5|. A tiny violation of isospin symmetry of parton distributions, largely 
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compatible with current data, would have a similar impact. Both these effects may have to be 
taken into account in the evaluation of the systematic error. 

• Generic corrections to the propagators or couplings of the SM gauge bosons can only 
produce a small fraction of the NuTeV anomaly, as shown in fig. In order to perform such a 
general analysis we have extracted the ‘oblique’ parameters and the SM gauge couplings directly 
from a fit of precision data, without imposing the SM predictions. We have assumed that the Z 
couplings are generation universal and respect SU(2)i, as in the SM. In principle, the NuTeV 
anomaly could be explained by new physics that only shifts the vZv couplings. However this 
situation is not realized by mixing the Z with extra vector bosons, nor by mixing the neutrino 
with extra fermions. 

• MSSM. Loop corrections in the MSSM have generally the wrong sign and are far too small to 
contribute significantly to the NuTeV observables. 

• Contact operators. Dimension six quark-quark-lepton-lepton operators can fit the NuTeV 

anomaly consistently with other data. The desired operators are neutral current, left-handed 
four fermion vertices of the form where q = {u,d}. The coefficient r] must 

be of order O.OI x 2^/2Gf, and the sign is fixed by requiring a negative interference with the 
SM. Effects of these operators should be seen at run II of Tevatron, unless they are generated by 
very weakly coupled light particles. If one restricts the analysis to SU( 2 ) 2 ,-invariant operators, 
only the operator [Qi7“(5i][T27aL2] can fit NuTeV. 

• Leptoquarks. singlet and triplet leptoquarks can induce these operators — but if 

the leptoquark masses are SU( 2 ) 2 , degenerate, either the sign is wrong or other unacceptable 
operators are also generated. A SU( 2 ) 2 , triplet leptoquark of spin one is a partial exception (see 
fig. 1^ at least from a purely phenomenological perspective. Non degenerate triplet leptoquarks 
could fit the NuTeV results, but squarks in i?-parity violating supersymmetry cannot. 

• Extra U(l) vector bosons. A Z' boson that does not mix with the Z boson can generate 
the NuTeV anomaly (see fig. ^), if its gauge group is B — 3L^, the minimal choice suggested 
by theoretical and experimental inputs. The Z' can be either heavy, 600 GeV ^Mz' ^5TeV, 
or light, 1 GeV <Mz' <10 GeV. The Z' that fits the NuTeV anomaly also increases the muon 
5-2 by 10 ® and (if light) gives a Z' burst to cosmic rays just above the GZK cutoff without 
requiring neutrino masses heavier that what suggested by oscillation data. 
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Note added The {v^) NuTeV beam contains a 1.7% contamination of (7'e)- A recent paper [^| 
suggests an explanation of the NuTeV anomaly assuming that 20% of oscillate into sterile 

neutrinos. The suggested oscillation parameters are not consistent with disappearance experiments 
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(unless Bugey and Chooz underestimated the theoretical error on their reactor fluxes). Furthermore, 
NuTeV not only predicts the and Ve fluxes through a MonteCarlo simulation, but also measures 
them directly (see pag. 26 of the transparencies in ^). The agreement between the two determinations, 
at the few % level, contradicts the oscillation interpretation. 


Note added after publication (June 2002) In a recent publication the NuTeV collaboration 
investigated the effect of a strange quark asymmetry and of isospin violation on their electroweak result. 
Specifically, they claim that the strange quark asymmetry is severely constrained by the dimuon data 
of Ref. |2^, and that the effect of strange quark asymmetry and isospin violation might be considerably 
diluted due to the fact that NuTeV does not measure directly the Paschos-Wolfenstein ratio Rpwj 
Eq. (I). _ 

The claim that the dimuon (i.e. vs —> /ic scattering) data p8|, 


strange asymmetry s 
s- = -(2.7 ± 1.3)10-3 


|67[| provides evidence against the 
+2 10-3 of the BPZ global fit ||^ (and in fact suggest a negative asymmetry 


|67[| ) appears dubious, for the following reasons: 


1. The parametrization of the strange and antistrange distributions assumed by NuTeV is un¬ 
physical, in that it violates the constraint that the proton and neutron carry no strangeness. 
Furthermore, due to its too small number of free parameters, it artificially relates the s/s asym¬ 
metry at a; < 0.5 (where NuTeV have significant data) to the one at x > 0.5 (where NuTeV have 
few events). 


2. NuTeV did not make a global fit allowing a strange asymmetry, but rather made a fit to their 
dimuon data, using a set of parton distributions based on pre-existing (now obsolete) fits, ob¬ 
tained neglecting NLO QCD corrections and optimized under the assumption s = s. This is 
especially worrisome due to the fact that the strange asymmetry found by NuTeV appears to 
depend very sensitively on the underlying set of parton distributions (see table I of |2^). 


I i3 


3. The dimuon data are considerably less sensitive to s than to s. In fact, the claim 

that a strange asymmetry at x > 0.5 is excluded at high confidence level should be restated as 
the statement that NuTeV rules out a total strangeness at x > 0.5 of the magnitude found by 
BPZ [p^]. However, what matters for the NuTeV anomaly is the strange asymmetry. 


The BPZ global fit |]^ is not subject to these drawbacks, but it did not include the recent dimuon 
data. The BPZ fit is characterized by a relatively large strange sea at large x, driven mainly by 
CDHSW data, which agrees well with positivity constraints derived from polarized DIS |^. Hence, 
the only conclusion that can be drawn by comparing these two analyses is that the size of the large-x 
strange sea suggested by CDHSW as analyzed by BPZ seems larger than allowed by NuTeV data p^ . 
The origin of this discrepancy, and its impact on the best-fit strange asymmetry, could only be assessed 
by performing a global NLO fit which includes all available data. Our statement in section 3 remains 
therefore unchanged: the impact of the data [^] on the strange asymmetry is unclear. 


NuTeV also comment on isospin violating parton distributions taking the model by Thomas 
et al. as reference. This model predicts a small effect on as a result of a subtle cancellation 
between high and low x regions. This conclusion is model-dependent. The fact remains that 0(1%) 
isospin violation effects could generate the NuTeV anomaly, without conflicting with any other existing 
data. 

Coming to the possible dilution of strangeness asymmetry or isospin violation due to experimental 
cuts, we should like to point out that, contrary to what stated in |67|, we did include charm threshold 
effects and some experimental cuts in our analysis. We found moderate dilution effects, as discussed 
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in section 3 of this paper. On the other hand, we cannot simulate the full NuTeV experimental set-up. 
However, if what NuTeV really measures differs from i?pw in n way which is significant at the required 
level of accuracy, then the cancellation of NLO effects that occurs in the Paschos-Wolfenstein relation 
cannot be taken for granted any longer. In particular, NuTeV claims to be less sensitive to Rp than 
is i?pw- In general any asymmetry between charged-current and neutral-current, or between v and 
V events spoils the cancellation of the NLO corrections in Eq. ( 0 ). Such asymmetries can also be 
induced by experimental cuts and by different v, D spectra. If any of these effects were significant, a 
NLO analysis would be required in order to obtain a reliable determination of at the desired level 
of accuracy. Only a full NLO analysis of the NuTeV data could settle this issue. 
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